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Motivation: 

The Launch Environment 



Significant resources have been spent to develop the materials and 
structure to withstand the harsh conditions of vehicle launches. 

The launch environment is highly complex in terms of geometric details 
and flow physics. 





Jet Impingement Problem: 

Pressure Environment 



Falcon Heavy 
pressure field 

launch 

platform 




flame deflector 


Motivated by CFD prediction 
requirements of launch pressure 
environments in the flame trench. 



damage after STS-124 launch 




Background: 

Jet Impingement Problem 

Supersonic jet impingement is an important problem for 
rocket launch and vertical take off and landing aircraft 

■ Complex instabilities when jet impinges on plate 

■ Tone generation for free jets is well documented 

(Tam (1990), Panda (1999), ...) 

■ Various studies on shock- wave oscillations 

(Henderson et al. (1993,2002,2005), Nakatogawa et al. (1971), Ginzburg et al. (1970), ...) 

■ Tone generation for certain nozzle pressure ratios and 
plate to nozzle distances 

(Krothapalli et al. (1999), Henderson et al. (2005), ...) 

■ Instability depends on plate location in the shock cell 
structure of corresponding free jet and standing-off shock 

(Henderson et al. (2005)) 

■ Feedback model: Downstream convected large scale coherent 
flow structures and upstream propagating pressure waves 
produce resonance 

■ Approximate formula for instability frequencies 

(Dauptain et al. (2012), Henderson and Powell (1993), Powell, Ho and Nosseir (1981), etc.) 



Exp.: Shadow photograph 

Henderson et al. (200 5) 

Sound 
wave 



LES: Density Grad 

Dauptain et al. (2012) 
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Approach: Post-Processing 

POD results in a decomposition of the flow field into a set of basis functions 
that capture most of the flow energy as defined by a user-defined norm with the 
least number of modes* 

i 

q(x,t ) ~ a^ n \t)x^ (x) 



n - 0 


*Rowley(200i), Freund and Colonius (2002) 
**Lumley(i967), Sirovich (1987), 

Chatterje (2004) 

***Nonomura et al (2011) 


Used snapshot method** in temporal domain 
Vector norm (energy) with | ^ 1 2 

q k =[p,u,v,w,T°-5] and co k =[i,o,o,o,o] 



^kRktfk ]dx 


k = 1 


DFT used to analyze time-spectral characteristics of pressure field 
amplitude ^9 spectral density peak frequency. 


4- 


Q 
"Kl 2 - 


0 - 




o 

y/D 



4- 


O 

2 - 


0 - 


l 




j 


0 

y/D 



HOSA results shown in paper (phase speed and coherence length) 



Approach: 

LAVA Solver 



Launch Ascent and Vehicle Aerodynamics 



■ Developed by authors at NASA Ames Research Center. 

■ Supports Cartesian AMR, block structured-curvilinear overset, and 
unstructured arbitrary polyhedral cells. 

■ Sharp interface immersed-boundary representation of geometry. 

■ Cell centered finite-volume and finite-difference formulation. 

■ Explicit and implicit time-integration. 

■ Multi-species capability with user-specified equations of state. 

■ Automatic volume gridding requiring only a surface triangulation. 


*Kiris et al.(2014), Brehm et al. (2014), Sozer et al. (2014), Housman et al. (2014), Moini-Yekta (2012,2013) 
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Approach: 

Investigate Level of Fidelity 
for Pressure Environment Simulations 

■ Slowly increase complexity of CFD model to examine assumptions: 

o Inviscid 

o Viscous (ILES, SA-DES), slip wall BCs* 
o Viscous (ILES, SA-DES), no-slip wall BCs 
-> What can these CFD approaches capture? 

■ Want fast turn-around times while maintaining high degree of accuracy 

o Reduce grid generation time 

-> immersed boundary methods essentially eliminate manual 
grid generation process 
o Reduce run time 

-> block structured mesh provides optimal memory layout 

*Dauptain et al., “Large Eddy Simulation of Supersonic Impinging Jets”, AIAAJ, 2012 
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Jet Impingement Problem: 

Computational Setup 

■ Collaboration with JAXA** 

■ Conditions based on experiments by 
Nakanishi et al. at the UT-Kashiwa 
hypersonic and high-temperature 
wind tunnel* 

■ M=i.8, air at T=300K (cold) 

■ Re=V e D/v=i.6xio 6 

■ 200-400 million grid points 

■ Nozzle-to-plate distance 5D 

■ Impingement angle a 

Flow conditions at nozzle exit and percent differences with perfectly expanded jet: 



Mach number 

pressure, p ref [Pa] 

exit velocity, V e [m/s] 

exit conditions 

1 .8045 

100,794 

488.14 

difference 

-0.25% 

-0.52% 

-0.24% 

*Nakanishi et al., “Acoustic characteristics of correctly-expanded supersonic jet impinging on an inclined at plate”, AJCPP2012-129 Jk 

**Tsutsumi et al. 2012, **Nonmura et al. 2010 & 2011 




Jet Impingement Problem: 

Unsteady Flow Field 



dilatation & |oo|D/V e «8.3 


pseudo-Schlieren 


i 







,J*A(p)/p re( df D/V e 


Jet Impingement Problem: 

Grid Convergence 



Spatial Resolution 



0.04 
0.03 : 

0.02 - 



Sub-Iterations 
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_i i 

20 40 

sub-iterations 


i 

60 


Point Probes 



# C = (0.55,3. 75)D 




# B = (5.0,0.0)D 


A = (4.3,-2.5)D 


■ Error convergence study w.r.t. spatial discretization, sub-iterations, time-step* 

■ Thorough verification study of LAVA in Moini-Yekta et al. **’*** 

■ Grid resolution sufficient until approximately St=3~4 

■ Much finer mesh than Dauptain et al. 


*Housman et al. “Space-Time Convergence Analysis of a Dual-Time Stepping Method for Simulating Ignition Overpressure Waves" ICCFD6, 2010^ 

**Moini-Yekta et al. “Verification and Validation Studies for the LAVA CFD Solver”, AIAA Paper, San Diego, June 23-27 th , 201 3 

***Kiris et al. (2014) 


D6, 2010. 


Jet Impingement Problem: 

Mean Flow 



Mean-Pressure Contours: 



■ Flow is mildly over-expanded 

■ Shock structure on the impingement surface 

■ Impingement pressure increases with angle 



Jet Impingement Problem: 

Mean Flow: 



Pressure Distribution in x=o cutting plane 

time-averaged approximated maximum 



■ Maximum in time-averaged pressure distribution analytically estimated (»3% error) 
(see for example Allgood for application to flame deflector) 

■ CFD needed to predict spatial pressure distribution and unsteadiness 

■ Maximum surface pressure is approximately 40-50% higher than time-averaged 
maximum surface pressure 
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Allgood, D. and Ahuja, V., “Computational Plume Modeling of Conceptual ARES Vehicle Stage Tests”, AIAA Paper 201 2 
Launcher/Transport Branch, “Development of Design Criteria for Saturn V Flame Deector", NASA-TM, 1965 


Jet Impingement Problem: 

Mean Flow: 



Pressure Distribution in x=o cutting plane 


time-averaged 
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■ Maximum in time-average | >timated (»3% error) 

(see for example Allgood for application to name ueneciorj 

■ CFD needed to predict spatial pressure distribution and unsteadiness 

■ Maximum surface pressure is approximately 40-50% higher than time-averaged 
maximum surface pressure 

Allgood, D. and Ahuja, V., “Computational Plume Modeling of Conceptual ARES Vehicle Stage Tests”, AIAA Paper 2012 
Launcher/Transport Branch, “Development of Design Criteria for Saturn V Flame Deector", NASA-TM, 1965 



Jet Impingement Problem: 

Mean Flow: 



Pressure Distribution in x=o cutting plane 

time-averaged approximated maximum 
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Allgood, D. and Ahuja, V., “Computational Plume Modeling of Conceptual ARES Vehicle Stage Tests”, AIAA Paper 201 2 
Launcher/Transport Branch, “Development of Design Criteria for Saturn V Flame Deector", NASA-TM, 1965 


Jet Impingement Problem: 

Inviscid vs. Viscous 


Max. surface pressure 



Streamlines and pressure rms 


mean pressure I I 

{■) / / / ] ! 
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ILES 


DES near the wall dampens pressure oscillations 

Very close agreement between viscous and inviscid results 


Jet Impingement Problem: 

Peak Frequencies 


Inviscid ILES 




DES 



■ General features of peak frequency distributions are very similar 

■ Several reasons for asymmetry in distribution: 

Noisy spectra, statistically not yet converged, (m=3)-mode, higher excitation 
from up-slope reflected waves 



Jet Impingement Problem: 

Peak Frequencies 


Inviscid ILES 




DES 



■ Several reasons for asymmetry in distribution: 

Noisy spectra, statistically not yet converged, (m=3)-mode, higher excitation 
from up-slope reflected waves 

■ General features of peak frequency distributions are very similar 



Jet Impingement Problem: 

Peak Amplitudes 

Inviscid ILES DES 

mm mm m m 



All three CFD models display similar distributions 

Reduced amplitudes for DES downstream of the impingement location 


Jet Impingement Problem: 

No-slip Wall BCs 



Streamlines 


Streamlines 


and mean pressure pressure rms 



■ No separation on impingement plane 


and mean pressure 

DES 

A i M i ' / •' 


pijimi 


pressure rms 




Note: curvilinear overset grid system 


■ Streamlines and pressure rms values are very similar 


■ Slightly more diffused solution due to coarser grid and metric terms 

■ Eddy viscosity from DES model seems to affect unsteadiness in shear layer 



Jet Impingement Problem: 





Jet Impingement Problem: 

Velocity Power Spectra (ILES) 



point probes 




■ Energy spectrum fills up from D to A (laminar-turbulent transition process) 

■ Points within red-shaded region show similar spectra 

■ Well resolved until approximately St= 3~4 

■ Frequencies correspond well to results in literature (Michaike (1984)) 
co@ c /U 0 «o. 1-0.2 for m=o and m=i (azimuthal modes) 




Jet Impingement Problem: 

Spectral Analysis: Amplitude 



■ Large amplitudes at impingement location - f low 
■ Broad spectrum in shear layer - f high 


0 0.05 o.i 

A 



Jet Impingement Problem: 

Spectral Analysis: 


Phase 



■ Provides overview of directionality of wave propagation 

o For acoustic waves as well as instability waves 

■ Shear layer waves close to nozzle lip appear to be unaffected by impingement 
angle for St=i.64, stronger effect for St=o.4i 





Jet Impingement Problem: 

Spectral Analysis: Phase 


y/D y/D y/D y/D 

Provides overview of directionality of wave propagation 
o For acoustic waves as well as instability waves 
Shear layer waves close to nozzle lip appear to be unaffected by impingement 
angle for St=i.64, stronger effect for St=o.4i 






Jet Impingement Problem: 

Proper Orthogonal Decomposition 



eigenvalues 


captured energy 
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Stronger coherence with increasing impingement angle 
~2o/ 50 modes capture 60/80% of energy (p) for 90° 


Jet Impingement Problem: 

POD modes 
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Jet Impingement Problem: 

POD modes 

45° 
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Jet Impingement Problem: 

POD modes 

45 ° 
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Jet Impingement Problem: 

Time Coefficients 
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St=fL D /V. 



Jet Impingement Problem: 

Time Coefficients 



<af 10 1 


mode 1 
mode 2 
mode 3 



St=fL D /V c 


St=fL D /V, 


■ Dominant frequencies, Sbo.36 and 0.33, for 45 0 and 90°, respectively 

■ Dominant frequency occurs for all three modes 



Jet Impingement Problem: 

Time Coefficients 



<af 10 1 



St=fL D /V e 


_%L L a^e. 


Dominant frequencies, St»o.36 and 0.33, for 45 0 and 90°, respectively 
Dominant frequency occurs for all three modes 




Jet Impingement Problem: 

Phase-Averaged IT- Velocity 



Phase-averaged at St *0.33 


Jet Impingement Problem: 

Spectral Analysis 



Downstream convected coherent structures interact with shock oscillations 
Generated acoustic waves excite shear layer 


Jet Impingement Problem: 

Spectral Analysis 




□OOUMp 

wave 


coherent 

structure 





Extract Fourier spectra along two lines 

(1) shear layer region 

(2) lxD away from nozzle lip 
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Jet Impingement Problem: 

Spectral Analysis 



Fourier spectrum 



ocoustic 

wave 


coherent 

ctructurc 



■ Extract Fourier spectra along two lines 

(1) shear layer region 

( 2 ) lxD away from nozzle lip 


Frequency (St»o.33) dominant for all 
vertical distances, z/D=i -5 
High frequencies in spectrum generated 
by linear and nonlinear mechanisms 
Low frequencies mainly due to shock 
oscillations 


ck^ 


Jet Impingement Problem: 

Spectral Analysis 



Fourier spectrum 



■ Extract Fourier spectra along two lines 
(2) lxD away from nozzle lip 


Dominant frequencies correspond to 
each other for (1) & (2) 


Jet Impingement Problem: 

Feedback Mechanism? 




Yes: 

■ Strong correlation between 
frequency spectra 


No: 

■ Unsteady coherent structures 
generate pressure waves 

■ Numerical noise 


Jet Impingement Problem: 

Feedback Mechanism? 



Pressure amplitude from DFT for St=o.4i 
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Yes: 

■ Strong correlation between 
frequency spectra 

■ Strong “Upstream” effect in 
pressure amplitudes values 


No: 0 0 05 ai 

■ Unsteady coherent structures 
generate pressure waves 

■ Numerical noise 




Jet Impingement Problem: 

Secondary Impingement 




Most fundamental studies focus on 
plate impingement without a 
horizontal plate 
Why is this important? 
o Most flame deflectors at KSC 
are shaped this way 
Flow physics are strongly affected 
by horizontal plate 
o Acoustics 

o Jet primary and secondary 
shock interactions 



Jet Impingement Problem: 

Secondary Impingement 



Mean pressure and streamlines 



rms pressure 
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y.D 

■ Strong unsteadiness at primary and secondary impingements 

■ Roughly 20 % larger peak pressure rms values 



Jet Impingement Problem: 

Secondary Impingement 


Mean pressure and streamlines 



rms pressure 
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No horizontal plate 



Strong unsteadiness at primary and secondary impingements 
Roughly 20% larger peak pressure rms values 


Jet Impingement Problem: 

Secondary Impingement 

Maximum pressure 



Mean pressure and streamlines 



no horizontal plate 



y/D 



y.D 

■ Strong unsteadiness at primary and secondary impingements 

■ Roughly 20% larger peak pressure rms values 



Jet Impingement Problem: 

Secondary Impingement (SI) 

Spectral Analysis: Phase 



Jet Impingement Problem: 

Secondary Impingement (SI) 

Spectral Analysis: Phase 



■ SI strongly affects acoustic near field 

■ SI seems to have affect on shear layer unsteadiness 
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Jet Impingement Problem: 

Secondary Impingement (SI) 





Spectral Analysis: Phase 


no horizontal plate 


St=o.82 



SI strongly affects acoustic near field 

SI seems to have affect on shear layer unsteadiness 




Jet Impingement Problem: 

Secondary Impingement 



Pressure amplitudes from DFT for St=o.4i 



■ Strong upstream effect of secondary impingement 

■ horizontal impingement plate is crucial for flame deflector analysis 



Jet Impingement Problem: 

Secondary Impingement 



Pressure amplitudes from DFT for St=o.4i 



■ Strong upstream effect of secondary impingement 

■ horizontal impingement plate is crucial for flame deflector analysis 



Jet Impingement Problem: 

Secondary Impingement 

Proper Orthogonal Decomposition (45 0 ) 
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Jet Impingement Problem: 

Secondary Impingement 

Proper Orthogonal Decomposition (30°) 
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Jet Impingement Problem: 

Computational Setup 




Cutting Planes 



Computational Setup: 

■ 50 billion grid points needed to obtain similar grid resolution as for model problem, 
150 million grid points used here ^23 grid points per nozzle diameter 

■ multi-species calculation, hot jet, inviscid, slip BCs 

■ 27 Merlin engines 




Jet Impingement Problem: 

Mean Pressure 
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Jet interaction 

Primary and secondary impingement pressure 
Step cannot be modeled with slip wall BCs 



Jet Impingement Problem: 

Pressure RMS Values 



High pressure rms values in Mach diamond structure 

Strong unsteadiness at primary and secondary impingement locations 


Jet Impingement Problem: 

Spectral Analysis 

peak frequency 
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flame deflector 


peak amplitude 
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■ Low frequency content at impingement locations 

■ High frequencies right above high pressure region 

■ Duct pressure oscillations due to containment 



Jet Impingement Problem: 

Proper Orthogonal Decomposition 



■ Interaction between jets and primary impingement 

■ 3 rd POD mode picks up secondary impingement 

■ Strong coherence suggest use of ROMs 

■ Coherence may be over-emphasized due to low grid 
resolution 
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Summary 

Computational Approach: 

■ slip wall BC valid assumptions for this problem 

■ Similar pressure signatures for inviscid, ILES, DES 

Jet Impingement Model Problem: 

■ Significant unsteadiness needs to be captured by CFD 

■ ILES captures unsteadiness until St«3-4 

■ Strong coherence in flow field: f low and f high 

■ POD provides insight about unsteadiness of shock oscillations 

■ 90°: Strong interaction of shock oscillations with unsteadiness of jet 

Feedback mechanism 

■ SI plays important role for pressure field on flame deflector 

■ SI effects unsteadiness of jet in near field 

Feedback mechanism 

Falcon Heavy: 

■ Similar flow physics as in impingement model problem 

■ Interaction between Mach diamond structure and shock oscillations 



Future Work 


Jet Impingement and Secondary Model Problems: 

■ Analyze transition process (azimuthal mode decomposition) 

■ Capture interaction between shock oscillation and jet flow by 
employing LST analysis 

■ 3D POD 

■ Use higher-order methods to capture higher-frequencies 

■ Reduced order model based on POD or disturbance flow formulation 


Falcon Heavy: 

■ Utilize ROMs to predict pressure field 


Thank you for your attention! 


STS-135 



Time = 0.007000 (s) 

exhaust gas mass fraction 



surface pressure contours 



